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Abstract 

Purpose As a consequence of the multi-functionality of land, 
the impact assessment of land use in Life Cycle Impact Assess¬ 
ment requires the modelling of several impact pathways cover¬ 
ing biodiversity and ecosystem services. To provide consistency 
amongst these separate impact pathways, general principles for 
their modelling are provided in this paper. These are refinements 
to the principles that have already been proposed in publications 
by the UNEP-SETAC Life Cycle Initiative. In particular, this 
paper addresses the calculation of land use interventions and 
land use impacts, the issue of impact reversibility, the spatial and 
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temporal distribution of such impacts and the assessment of 
absolute or relative ecosystem quality changes. Based on this, 
we propose a guideline to build methods for land use impact 
assessment in Life Cycle Assessment (LCA). 

Results Recommendations are given for the development of 
new characterization models and for which a series of key 
elements should explicitly be stated, such as the modelled land 
use impact pathways, the land use/cover typology covered, the 
level of biogeographical differentiation used for the character¬ 
ization factors, the reference land use situation used and if 
relative or absolute quality changes are used to calculate land 
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use impacts. Moreover, for an application of the characterisa¬ 
tion factors (CFs) in an LCA study, data collection should be 
transparent with respect to the data input required from the 
land use inventory and the regeneration times. Indications on 
how generic CFs can be used for the background system as 
well as how spatial-based CFs can be calculated for the 
foreground system in a specific LCA study and how land 
use change is to be allocated should be detailed. Finally, it 
becomes necessary to justify the modelling period for which 
land use impacts of land transformation and occupation are 
calculated and how uncertainty is accounted for. 

Discussion The presented guideline is based on a number of 
assumptions: Discrete land use types are sufficient for an 
assessment of land use impacts; ecosystem quality remains 
constant over time of occupation; time and area of occupation 
are substitutable; transformation time is negligible; regenera¬ 
tion is linear and independent from land use history and 
landscape configuration; biodiversity and multiple ecosystem 
services are independent; the ecological impact is linearly 
increasing with the intervention; and there is no interaction 
between land use and other drivers such as climate change. 
These assumptions might influence the results of land use Life 
Cycle Impact Assessment and need to be critically reflected. 
Conclusions and recommendations In this and the other 
papers of the special issue, we presented the principles and 
recommendations for the calculation of land use impacts on 
biodiversity and ecosystem services on a global scale. In the 
framework of LCA, they are mainly used for the assessment 
of land use impacts in the background system. The main areas 
for further development are the link to regional ecological 
models running in the foreground system, relative weighting 
of the ecosystem services midpoints and indirect land use. 

Keywords Biodiversity • Ecological functions • Ecosystem 
services • Life cycle assessment • Land use 

1 Introduction 

Life Cycle Assessment (LCA) is a tool to support decision¬ 
making widely used to assess the potential environmental 
impacts of a given product/service at each step of its life cycle. 
Impacts at the global and regional scales are accounted for, 
such as climate change, eutrophication, acidification, toxico¬ 
logical effects and abiotic resource use. Early attempts to 
integrate direct and indirect land use in LCA (Lindeijer et al. 
2002; Mila i Canals 2007a; Reinhard and Zah 2009) and its 
impact on biodiversity (e.g. Koellner and Scholz 2007; 
Koellner and Scholz 2008; Michelsen 2008; Schmidt 2008; 
Geyer et al. 2010a) and on ecological functions (e.g. Lindeijer 
et al. 2002; Mila i Canals L 2003; Mila i Canals L 2003; Mila i 
Canals 2007a; Beck et al. 2010; Saad et al. 2011) have been 
made, but these are still not fully operational or widely applied. 


Although functional species diversity is an important factor in 
the cause-effect chain from land use to ecosystem functioning 
and services (Balvanera et al. 2006; Flynn et al. 2009), func¬ 
tional aspects of biodiversity are not yet considered (see the 
review from Curran et al. 2011). Moreover, only preliminary 
studies explicitly discuss ecosystem services in LCA (Maes et 
al. 2009; Zhang et al. 2010a, b; Bare 2011) and this concept 
generally refers to life support functions (de Groot 1992; de 
Groot et al. 2002) introduced early into LCA (e.g. Udo de Haes 
et al. 2002; Anton et al. 2007; Mila i Canals 2007b). 

The problem with all those attempts to integrate land use 
impacts on biodiversity, ecological functions and ecosystem 
services into LCA is their limited focus on one geographical 
scope such as case study regions of some square kilometres 
or larger specific biogeographical regions like Canada or 
Central Europe. However, the strength of LCA is to provide 
a life cycle perspective (Mila i Canals et al. 2007) and thus 
requires methods that are able to account for an assessment 
of land use impacts related to a variety of land use types and 
locations. Even simple products like milk imply globally 
distributed land use through e.g. the supply chain of the 
concentrate feed (Cederberg and Mattsson 2000). Also the 
increasing demand for studies on biofuels and their environ¬ 
mental footprint shows that a globally applicable Life Cycle 
Impact Assessment (LCIA) method is needed to assess land 
use impacts on biodiversity and ecosystem services. 

In summary, land use impact assessment in LCA has 
moved significantly since the early works of SETAC 
(Lindeijer et al. 2002) and the first phase of the UNEP- 
SETAC Life Cycle Initiative (Mila i Canals et al. Mila i Canals 
2007a), with increasing attempts to integrate land use impacts 
in LCA. Now, more consistency is needed in order to take 
stock of what has been suggested to date and provide align¬ 
ment in the future modelling of different impact pathways. 

This paper builds on the “Key Elements in a Framework 
for Land Use Impact Assessment Within LCA” that was 
developed in the context of the UNEP-SETAC Life Cycle 
Initiative (Mila i Canals 2007a). The aim of the current 
paper is to move beyond the description of key elements 
by suggesting specific guidelines for a comprehensive and 
consistent impact assessment encompassing all pathways 
that originate from land use and damages on biodiversity 
and ecosystem services. Such guidelines are applied and 
exemplified throughout this special issue, first by develop¬ 
ing a globally consistent Life Cycle Inventory (LCI) classi¬ 
fication system (Koellner et al. 2013) then, by proposing 
new methods for the different impact pathways affecting 
biodiversity (de Baan et al. 2013; Souza et al. 2013) and 
ecosystem services (Muller-Wenk and Brandao 2010; Beck 
et al. 2010; Brandao and Mila i Canals L 2013; Saad et al. 
2013) and finally, by illustrating the previous new methods 
and recommendations in a case study on margarine produc¬ 
tion (Mila i Canals L et al. 2013). 
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2 Principles of a globally applicable land use impact 
assessment method 

2.1 Land use interventions and their impacts on ecosystem 
quality 

Two types of land use interventions are usually considered 
in life cycle inventories and impact assessments; land trans¬ 
formation and land occupation (Lindeijer et al. 2002; Mila i 
Canals 2007b). During land transformation (also called land 
use change, LUC), the properties of a piece of land are 
modified to make it suitable for an intended use such as 
deforesting or draining land to establish arable fields. The 
phase of transformation is relatively short, and the temporal 
dimension is neglected. During land occupation, land is 
used in the intended productive way (e.g. arable field) and 
the properties of a piece of land are maintained (e.g. the 
regrowth of forest is avoided on an arable field). 

These land use interventions have an impact on ecosys¬ 
tem quality Q over a certain period of time, where Q can be 
defined as the capability of an ecosystem (or a mix of 
ecosystems at the landscape scale) to sustain biodiversity 
and to deliver services to the human society. This refers to 
the area of protection natural environment , which provides 
the “intrinsic value of nature (ecosystems, species) and the 
economic value of life support functions” to the human 
society (Udo de Haes et al. 1999a, b). Land use impacts 
and damages to ecosystem quality may be measured with 
different indicators expressing the intrinsic value of biodi¬ 
versity and natural landscapes or the functional value of 
ecosystems in terms of their goods (i.e. natural resources 
like timber or food) and services (i.e. life support functions 
like climate regulation or erosion regulation) (sensu Mila i 
Canals 2007a). These impacts result from both land occu¬ 
pation (because ecosystem quality is kept at a different level 
than would naturally/otherwise be present) and land trans¬ 
formation (because the characteristics of ecosystems are 
changed on purpose). 

Fig. 1 Simplified illustration of Ecosystem quality Q 
transformation impact (77) and 
occupation impact (07) for 
three land use types with 
different regeneration rates 

( J'LUl, reg? Cc/2, reg 3nd A US, reg)- 

For simplicity, the area A of 
occupation or transformation, 
which would embrace the third 
dimension, is not shown in the 
graph, but in the equations 


If no occupation process would follow a land transfor¬ 
mation, the forces of nature would gradually change the 
ecosystem quality towards its initial quality (although the 
original ecosystem quality might not be reached, see Section 
2.2). The impact of land transformation (77) is calculated as 
the integral of the difference in ecosystem quality between 
the land use situation and a suitable reference (A Q) over 
time, multiplied by the transformed area (A). As the tempo¬ 
ral dynamics of ecosystem quality are mostly unknown, a 
linear trajectory of ecosystem regeneration is assumed 
(Fig. 1). Transformation impacts are then calculated using 
Eq. (1), whereby t reg is the time required for full regenera¬ 
tion of ecosystem quality. 

77 = 0.5 x A Q x t reg x A (1) 

The inventory flow records the area A transformed and 
the characterization factor of transformation CF trans is given 
in Eq. (2): 

CF trans — 0.5 a AQ x t reg (2) 

Accordingly, the impact of land occupation ( OI) is cal¬ 
culated as the integral of A Q over time multiplied by the 
occupied area A (Eq. (3)). By assuming that A Q is constant 
during the occupation phase, occupation impacts are calcu¬ 
lated following Eq. (3) (with 7 being the time that a piece of 
land is occupied). 

OI = AQxT xA (3) 

Here, the inventory flow is given as A occ x T occ and the 
characterization factor of occupation CF occ is given in Eq. 
(4): 

CF occ = AQ (4) 

Figure 1 illustrates three examples of typical land use 
interventions and their associated impacts on ecosystem 
quality. For simplicity, the area A of occupation or transfor¬ 
mation, which would embrace the third dimension, is not 
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shown and the measured impact indicator discussed refers to 
impacts on biodiversity. However, the same reasoning ap¬ 
plies for other land use impact pathways. First, at time t\, the 
land is transformed from a reference (e.g. broadleaf lowland 
forest) to land use type LU1 (e.g. species-rich dry meadow), 
indicating a higher ecosystem quality (if measured as biodi¬ 
versity) than the reference. The transformation impact is 
given as the difference in ecosystem quality ( Qrej^Quji ) 
multiplied by the time it would take after abandoning LU1 
to restore the reference (e.g. the time required for a lowland 
broadleaf forest to naturally establish on a dry meadow). 
Both the transformation (area I in Fig. 1) and occupation 
impact (II) result in negative values, which denote a benefit 
for ecosystem quality. In the second situation, at time f 3 , 
land is transformed from a reference (e.g. broadleaf lowland 
forest) to land use type LU2 (e.g. intensive pasture). Here, 
the land transformation (III) and occupation (IV) show 
damaging impacts on ecosystem quality (i.e. positive values 
for 77 and OI). In the third situation (time f 4 ), land is 
transformed from LU2 (e.g. intensive pasture) to LU3 (e.g. 
intensive arable crops). The occupation impact (VI) is cal¬ 
culated as in the two previous examples. The transformation 
impact (area V) can be calculated by subtracting the impact 
of transforming a land from the reference to LU2 ( TI re f _ >LU2 ) 
from the impacts of transforming land from the reference to 
LU3 ( TI re f _ >LU3 ; see Eq. V in Fig. 1). Also, it is important to 
note that the regeneration times of ecosystems are dependent 
on the land use type (t LU1>reg ±t LU1>reg ±t LU1>reg ). Details on 
the separation of transformation into two separate flows 
“transformation from” and “transformation to” are given in 
the Electronic Supplementary Material (Figs. SI, S2 and 
S3). 

2.2 Reversibility of impacts from land use and permanent 
impacts 

In the above reasoning, it is assumed that after a certain 
regeneration time, the ecosystem quality of the reference 
situation could be re-established. Of course, an ecosystem 


that has been changed by human activities or by natural 
forces will never be exactly the same. Ecosystem quality is a 
result of the interaction between life (ecosystems) with the 
abiotic environment, and life is, in the strict sense of the 
word, not reversible. However, the forces of nature are able 
to veil the traces of most human activities if these forces are 
free to act during decades, centuries or millennia: aban¬ 
doned human settlements or structures are sooner or later 
wiped out and reconverted to a quasi-natural land cover 
which depends on the biogeographical conditions of the 
location. 

In the context of land use, it is proposed to consider 
impacts generally as reversible in the broad sense of revers¬ 
ibility. This means that abandoned land spontaneously de¬ 
velops towards a site-dependent potential natural vegetation 
(PNV) if the absence of human action continues during a 
sufficient length of time (regeneration time, also called 
relaxation time). At the end, an abandoned area can be 
considered as roughly equivalent, although not identical, to 
its pre-impact state. However, there are situations where the 
regeneration time, according to current knowledge, will 
exceed the modelling horizons of usual LCA studies, or 
even will exceed any finite number of years: a high salinity 
area in very dry climate could be barren for an indefinite 
time period. Such impacts are called permanent impacts. 
Permanent impacts can be expressed by multiplying the 
difference A Q between the initial reference ( Q re j ) and a 
new established steady state ( Q re f 2 ) by the area of transfor¬ 
mation and, alternatively, a certain modelling time (see 
equations in Fig. 2). Permanent impacts could also be quan¬ 
tified without choosing an (arbitrary) time horizon. Howev¬ 
er, it is recommended to multiply with a certain modelling 
time to get a standard unit of A Qx time x area for all land use 
impacts, which allows to easily aggregate permanent, trans¬ 
formation and occupation impacts in the same units. 

Nevertheless, such a pragmatic aggregation may have 
several issues regarding interpretation and effect on results. 
On the one hand, permanent impacts represent diminishing 
options for future development of a piece of land. On the 


Fig. 2 Calculation of 
permanent impacts caused by 
land use change. For simplicity 
again, the area A of occupation 
or transformation, which would 
embrace the third dimension, is 
not shown in the graph, but in 
the equations 


Ecosystem quality Q 



Calculation of impacts: 

Permanent impact: 

P'ref^Lui = (Gref- Q r ef2 ) *t M * A (equals area lb) or 
PI ref -, LU1 = (Qref~ G ref2 ) * A (delta between quality states) 


Occupation and transformation impact: 

0 ~P'ref >LU1 = [0-5 *(Q ref2 ~ Q LU1 ) * (t 3 -t 2 ) 

+ (Q ref - Q re f 2 ) *t M ] * A (equals areas la plus lb) 

ll)OI LU1 = (Q ref2 - Q LU1 ) * (t 2 -t 2 ) * A 

IP) TI L U1^U2= [0-5* (Q ref2 — Q LU2 ) * (t 5 -t4) 

- PI ref >LU1 ] A 

IV) OI LU2 = (Q ref2 - Q LU2 ) * (t 4 -t 2 ) * A 
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other hand, occupation and transformation impacts rather 
describe actual, temporary impacts occurring during 
occupation/regeneration phase. For this reason, it can be 
argued that aggregation of temporary and permanent im¬ 
pacts is equivalent to aggregation of different impact cate¬ 
gories. This implies weighting them against each other, but 
the value judgement, expressed by the time horizon over 
which permanent impacts are considered, should be made 
explicit as a value choice. As an overall consensus has not 
been reached among this workgroup, it is highly recommended 
to perform sensitivity analyses and discuss the influence of the 
modelling choices on the conclusion of the study, i.e. aggre¬ 
gating or not permanent impacts into a single result and, in case 
of aggregation, the choice of timeframe. 

Figure 2 illustrates how permanent, transformation and 
occupation impacts are calculated if full recovery of ecosys¬ 
tem quality is not possible within the modelling time. Ini¬ 
tially (time * 1 = 0 ), land is transformed from a reference (e.g. 
tropical rain forest) to land use type LU1 (e.g. pasture), 
having a permanent (damaging) effect on ecosystem quality. 
After a regeneration period of t 3 -t 2 years, a new steady state 
Qref 2 (e.g. old growth secondary forest) is reached. Trans¬ 
formation and occupation impacts of LU1 as well as the 
subsequent LU2 (e.g. arable land) are then calculated based 
on A Q=Q ref 2 ~Qlui/lu2 (see Eqs. I-IV in Fig. 2). The 
permanent impact can either be expressed by the A between 
Q re f and Q re f 2 , represented by the curled bracket, or by the 
area lb (see Fig. 2), representing the maintenance of an 
ecosystem quality below the initial reference (Q re f) for the 
modelling period defined (e.g. t M = 500 years). Permanent 
(and transformation) impacts would then be allocated to the 
products obtained from land during the first 20 years (see 
Section 3.2.4). 

Based on this approach, impacts occurring after the 
modelling period (i.e. Q re f~Qrep after year 500 in Fig. 2) 
are not accounted for. However, it is important to note that 
such a modelling decision is similar to the modelling 
choices adopted by the modelling period of global warming 
potential (see Section 3.3.1). 

2.3 Impact proportional to relative or absolute quality 
change? 

As mentioned before, the magnitude of permanent, transfor¬ 
mation and occupation impacts is dependent on the differ¬ 
ence in ecosystem quality (A Q) between a reference and a 
land use situation. A Q can be calculated based on absolute 
or relative differences, and which of both is preferable is 
subject of an ongoing debate. 

For the land use impact on carbon storage, the absolute 
A Q is clearly preferable as the goal is to identify the amount 
(in tonnes) of sequestered C released into the air in form of 
C0 2 due to the land transformation. This contributes to 


impact on climate change (and consequently on human health 
and non-human life), which is controlled by absolute flows in 
C or C0 2 tonnes. Measuring A Q in relative terms would not 
make sense, as releasing 50 % of the carbon stored in a carbon- 
rich vegetation has clearly not the same impact as releasing 
50 % of the carbon stored in a carbon-poor vegetation. The 
same holds true for other ecosystem services, as mostly the 
absolute amount of service provided is of importance. 

Less clear is the case of the biodiversity impact. If a 
factory building is erected on species-poor land (e.g. with 
10 species), or alternatively on a species-rich land (e.g. with 
100 species), the relative biodiversity A Q would be 100 % 
for both alternatives, whilst the impact on species-rich land 
is ten times higher than on the species-poor land in absolute 
terms. Selecting relative or absolute impacts is finally based 
on a value choice: choosing relative impacts gives equal 
weight to ecosystems, whereas choosing absolute impacts 
gives equal weight to species. 

In conclusion, it is recommended to calculate absolute 
impacts for ecological functions and ecosystem services on 
a global scale. When compared to relative impacts, such 
results are easier to interpret and allow coordination be¬ 
tween different impact pathways. However, both absolute 
and relative impacts are deemed appropriate and are a matter 
of value choices for biodiversity. To allow better compari¬ 
son of results across impact pathways, assessing biodiversi¬ 
ty in absolute terms would be strongly advised. However, 
the challenge remains as current heterogeneity of available 
empirical data on a global level do not necessarily allow for 
an absolute impact assessment. 

3 Results: guideline of a global land use impact 
assessment method 

Building on the recommendations by (Mila i Canals et al. Mila 
i Canals 2007a), this paper provides guiding principles for the 
development of impact assessment methods for biodiversity 
and ecosystem services. The guideline suggests options with 
respect to the creation of a spatial model, the inventory data 
collection as well as the land use impact calculation (Fig. 3). 
Based on the work of Saad et al. (2011), Fig. 3 highlights a 
series of key elements to be accounted for when developing a 
land use impact assessment method and thus proposing new 
characterisation factors (CFs). 

In an effort to provide a transparent and comprehensive 
approach for the creation of the model, it is suggested to 
explicitly state and define: which impact pathways with 
respect to biodiversity and ecosystem services are modelled 
(la), which land use/cover typology (lb) as well as the 
biogeographical differentiation level (lc) are used for the 
development of CFs and, in addition to the reference situation 
(Id), whether relative or absolute quality changes (le) are 
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Fig. 3 Elements of the UNEP- 
SETAC guideline to build a 
land use impact assessment for 
biodiversity and ecosystem 
services (adapted from Saad et 
al. 2011) 





used for the calculation of land use impacts. In practice, when 
applying CFs in an LCA study, the data collection should be 
transparent with respect to the data input required from the 
land use inventory (2a) and the regeneration time defined. It is 
also recommended to indicate how generic CFs can be used 
for the background system as well as how spatial-based CFs 
can be calculated for the foreground system if needed in a 
specific LCA study (2c) and how land use change is allocated 
to functional units (2d). With regards to the land use impact 
calculation, it becomes necessary to justify the modelling 
period for which the impacts of land transformation and 
occupation are calculated (3 a) and finally how the uncertainty 
of the impact assessment is assessed (3b). 

The following sections provide guidance on how to apply 
such principles, and examples of their application are pro¬ 
posed in the methodological papers of this special issue. 

3.1 Creation of spatial model 

3.1.1 Modelled impact pathways 

Land use intentionally and unintentionally influences the 
biodiversity as well as the structure and functions of ecosys¬ 
tems, causing damages to the areas of protection as defined 
in Jolliet et al. (2004). Figure 4 shows the cause-effect chain 
linking the land use with impacts on the four areas of 
protection. In order to enhance global relevancy of LCIA 
research, it is essential to link LCA activities to other ongoing 
research on ecological impact assessment. For this reason, a 
structure of the LCIA, which is in accordance with the 
globally acknowledged typology of ecosystem services 
of the Millennium Ecosystem Assessment (MA 2005) is 


proposed. Therefore, two main impact pathways are distin¬ 
guished: biodiversity damage potential and ecosystem services 
damage potential (Table 1). For each of the impact pathways, 
Table S2 (Electronic Supplementary Material) holds Excel files 
with the generic CFs. The calculation of those CFs is described 
in detail for biodiversity (de Baan et al. 2013; Souza et al. 2013) 
and for the ecosystem services (Muller-Wenk and Brandao 
2010; Brandao and Mila i Canals 2013; Saad et al. 2013). 

The biodiversity damage potential includes the protection 
of global species diversity and also the functional diversity 
of species in ecosystems. The ecosystem services damage 
potential is structured according to the classification 
suggested by the Millennium Ecosystem Assessment into 
services provided by ecosystems. The focus lays on the 
development of generic characterization factors for the fol¬ 
lowing impacts linked to land use: the impact on the poten¬ 
tial of the ecosystem to produce biomass (Biotic Production 
Potential); the impact on climate by influencing the carbon 
sequestration in the top soil and land cover (Climate Regu¬ 
lation Potential); the impacts on water quantity and quality 
(freshwater regulation potential and water purification po¬ 
tential); and the impacts on soil quantity and quality (erosion 
regulation potential). This is an initial list of recommended 
impact categories linked to land use for which characteriza¬ 
tion models were available or under development. Nonethe¬ 
less, the framework remains open to include additional 
midpoint oriented impact indicators if needed, such as the 
microbial activity indicating soil fertility. This would only 
be necessary if the new indicators show different results 
than those already provided by the current suggested list. 
In a comprehensive approach, new characterization methods 
for land use LCIA should specify the impact pathways 
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Fig. 4 Cause-effect chain for land use impacts on biodiversity and ecosystem services (adapted from Lindeijer et al. 2002). For more explanation 
on the specific impact pathways, see Table 1 and the papers in this special issue 


modelled within the framework. Future studies exploring 
the links between different ecosystem services, and 
between ecosystem services and biodiversity are strongly 
recommended. 

3.1.2 Land use and cover typology 

For a globally applicable land use LCIA, it is desirable that land 
cover classes and land use types related to them are determined 
in a consistent and generally accepted way for all continents, 
and that available land cover data are recorded according to 
such a classification. A hierarchical approach of land use clas¬ 
ses in LCA is given in this special issue (Koellner et al. 2013). It 
consists of four levels of detail, ranging from very general 
global land cover classes, to more refined categories and very 
specific categories indicating land use intensities. 

New LCIA methods should be able to differentiate between 
all major land use types generally accounted for in LCA and 
specify the level of differentiation achieved. Thus, the classifi¬ 
cation suggested by Koellner et al. (2013) can be used as a 
starting point. It is foreseeable that the same level is not achieved 


for all major land use types (e.g. a method may differentiate 
agricultural land up to level 3, while other land uses such as 
forestry or artificial areas are only specified at level 1). 

3.1.3 Bio geo graphical differentiation in land use impact 
assessment 

Biogeographical differentiation in land use LCIA is a major 
aspect to be considered since the same type of land use may 
trigger different impacts at different locations of the globe 
(Jolliet et al. 2004). Koellner et al. (2013) in this special 
issue suggest a hierarchical regionalisation on a global scale 
based on the classification system of biomes and ecoregions 
(Olson et al. 2001). Similarly to the case of land use typology, 
new LCIA methods for land use impacts should specify the 
level of biogeographical differentiation achieved. 

3.1.4 Reference situation 

As explained in Section 2.1 land use impacts as assessed in 
LCA are proportional to the difference in quality between 
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Table 1 Operationalisation of characterization factors in Life Cycle Impact Assessment for land use 
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the studied system and a reference situation. Three main 
options are proposed to describe the reference: 

• Option 1 is the concept of PNV, “which describes the 
expected state of mature vegetation in the absence of 
human intervention.” (Chiarucci et al. 2010). 

• Option 2 is to take as the reference states the (quasi-) 
natural land cover in each biome/ecoregion, i.e. the 
natural mix of forests, wetlands, shrubland, grassland, 
bare area, snow and ice, lakes and rivers. 

• Option 3 is the current mix of land uses as reference as 
proposed for Europe (Koellner and Scholz 2008). 

Options 1 and 2 are probably close enough to each other for 
the purposes of LCA. On the one hand, global data on the 
properties of potential “natural” land for many biogeographical 
regions are available in more or less satisfactory quality, so that it 
is globally practicable to assess anthropogenic land use impacts 
against the potential “natural” land cover defined as a reference. 
On the other hand, the current mix of land uses is a moving 
yardstick and would not be very practical in terms of definition. 

In conclusion, and in accordance to Mila i Canals et al. 
(2007a), it is recommended to use the (quasi-)natural land 
cover predominant in global biomes and ecoregions as a ref¬ 
erence when assessing land use impact on a global scale. 
Nevertheless, defining a reference situation is an area for 
further exploration recognised as a value choice, as using 
options 1 or 2 vs. option 3 supports different types of decisions 
(see Mila i Canals L et al. 2013). 

3.1.5 Absolute or relative land use impacts 

For each land use LCIA, it should be made explicit if land use 
impacts are assessed as absolute or relative changes to a 
reference situation. However, based on the above reasoning 
(see Section 2.3), it is generally recommended to calculate 
absolute land use impacts. 

3.2 Data collection 

3.2.1 Land use inventory data 

Land use inventory flow for land occupation and land trans- 
formation are measured in m years and in m , respectively. 
Although the level of differentiation (land use typology and 
spatial scale) may vary between methods (see Table SI in 
the Electronic Supplementary Material), both the land use 
type and biogeographical information should be registered 
based on the proposal on global land use inventories 
presented in this special issue (Koellner et al. 2013). The 
level of detail in the land use type classification and the 
ecosystem type refinement depends on the scope of the LCA 
study and on whether it occurs in the foreground or the 


background system; see e.g. the discussions in the case 
study presented in this special issue (Mila i Canals L et al. 
2013). If land use types are completely unknown in the 
product system, assuming a mix of intensive land use com¬ 
posed of forest, pasture/meadow, (non-) irrigated arable, 
permanent crops and artificial areas based on the propor¬ 
tions in the Global Land Cover 2000 is recommended. To 
assume intensive land use should stimulate the search for 
specific information for unknown land use types. 

3.2.2 Regeneration time 

Assuming that land use impacts are reversible (reversibility 
in the broad sense), the key question is, how many years 
after abandonment are needed until man-made traces in 
vegetation and soil have essentially disappeared? The length 
of this regeneration time, essential for the calculation of 
transformation impacts (see Section 2.1), depends mainly 
on the following factors: 

• The impact pathway: on a given plot, it may take more 
time to regenerate biodiversity than to regenerate the 
biotic production capacity 

• The type of land transformation: it takes more time to 
regenerate a forest, if the transformation resulted in 
sealed urban land, in comparison to a transformation 
resulting in cropland 

• The biogeographical conditions of the location: regen¬ 
eration is generally faster in warm and humid climate 
than in cold or dry climate. 

As a consequence, each of the impact pathways requires 
a set of regeneration times per type of transformation and 
per climatic region. However, as current land use intensity 
in many parts of the world is a historically new phenome¬ 
non, knowledge on regeneration times of ecosystems within 
intensely used landscapes is limited or hardly available. 

For impact pathways that are mostly dominated by 
soil carbon, it is recommended to use the regeneration 
time values for carbon storage in vegetation and soil 
provided by Muller-Wenk and Brandao (2010) for the 
impact pathway carbon sequestration potential (CSP). 
However, one should acknowledge that regeneration 
times for other impact categories derived from CSP data 
in this way are highly uncertain and strongly influence 
the transformation results. 

With respect to the biodiversity damage potential, 
Koellner and Scholz (2007) give rough estimates for 
regeneration times based on Bastian and Schreiber 
(1999). Those data indicate that mature peat bogs and 
old growth forests take up to 1,000 years for complete 
regeneration in the temperate zone. This regeneration of 
forests might be speeded up in tropical climate to 100- 
300 years if species richness is taken as an indicator, 
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but roughly the same when endemism is considered as 
the benchmark (Liebsch et al. 2008). 

Generally, regeneration time increases towards poles and 
with altitude. It also depends on the impact type as well as 
the ecosystem function or taxonomic group considered 
(Jones and Schmitz 2009). To take this into consideration, 
data based on altitude and latitude made for regeneration of 
biomass after land has been cleared, assuming no soil deg¬ 
radation are provided in (van Dobben in Lindeijer et al. 
1998, Annex 1.23). These numbers can be used as proxies 
to estimate regeneration of cleared land in different biomes 
or ecoregions, by calculating an area weighted average per 
region (depending on the share of total land within different 
altitudinal and latitudinal classes). However, further as¬ 
sumptions are needed to estimate regeneration time of dif¬ 
ferent land use types within one biome or ecoregion. For 
example, a constant regeneration rate r (decrease of A Q per 
time) could be assumed across land use types, which can be 
calculated as: 


r AQ/ t reg ( Q re f Q cleared ) / Leg, cleared 


( 5 ) 


where Q re f is the ecosystem quality of the reference, Qdeared 
that of a land from which vegetation cover was removed 
completely but no soil degradation took place, and t reg> cleared 
the time needed to regenerate an ecosystem after it has been 
cleared. However, such assumptions clearly do not reflect 
reality and should be used with caution. Since results of 
overall impact calculation for land transformation are very 
sensitive to regeneration time (see Schmidt 2008) and this 
factor is highly uncertain, we recommend providing uncer¬ 
tainty estimates by LCIA developers and using sensitivity 
analysis when applying CFs (e.g. using low and high esti¬ 
mates of regeneration times). 


3.2.4 Allocation of land transformation impacts 

Impacts from land transformation have to be allocated to 
output (functional units) arising from the new land use. This 
allocation or amortisation of land transformation belongs to 
the inventory, whereas calculating the magnitude of the land 
transformation impacts (integration of A Q between studied 
land use and reference over time, per area of land transfor¬ 
mation) lies in the LCIA phase. 

One option is to take a fixed amortisation period for 
all land use types. Twenty years would be consistent 
with IPCC suggestion for soil organic carbon emissions 
(IPCC 1996; Flynn et al. 2012). Allocating the transfor¬ 
mation impact on the production output of 20 years 
represents a good compromise between allocating them 
all to the first year (and thus quickly loosing sight of 
the effects of LUC) and using a long allocation period 
(which could lead to a quasi-elimination of transforma¬ 
tion impacts in the LCA results). An alternative is to 
assume different amortisation periods depending on the 
land use type. This was proposed in Ecoinvent v 2.0 
(Frischknecht and Jungbluth 2007, Table 5.5) based on 
default use periods. Alternatively, a linear depreciation 
along the regeneration pathway could be applied. For 
example, Fig. 1 shows that transformation impact TI re f . 
>lu2 would be 0 when f is reached. After that, the 
impact is 100 % related to the occupation. This means 
the regeneration time for a specific land transformation 
equals the appreciation period. 

However, the actual number of years is arbitrary in any case. 
In the absence of a clear, scientifically robust alternative, we 
suggest using 20 years as an allocation period, as this is in line 
with standards and regulations for land use-derived greenhouse 
gas emissions allocation (IPCC 1996; BSI 2008) but also 
recommend sensitivity analysis. 


3.2.3 Generic characterization factors versus case 
dependent calculation of impacts 

The use of generic characterization factors in LCA is a 
simplification in order to allow the calculation of land 
use impacts for the background system. As land use is 
often globally distributed, but the impacts on biodiver¬ 
sity and ecosystem services are clearly depending on the 
type and location of the land use, we recommend using 
differentiated CFs and link them consistently to differ¬ 
entiated land use inventories as proposed in this special 
issue. However, for the foreground system specific, CFs 
can be calculated and should be applied in LCA as 
shown by the LANCA tool for ecological functions 
(Beck et al. 2010) or by Geyer et al. (2010a, b) for 
biodiversity. This should enhance the validity of the 
impact assessment and should also reduce uncertainties. 


3.3 Land use impacts calculation 

3.3.1 Modelling period for land use impacts 

The basic principles for the calculation of land use impacts 
from occupation and transformation are shown above (Section 
2.1). Occupation impacts are modelled over the occupation 
time and thus no consideration of modelling period is neces¬ 
sary. But for transformation impacts, this calculation requires 
defining the modelling period: “the time lag between the land 
use intervention and the impact may be large; thus land use 
impacts should be calculated over a reasonable time period 
after the actual land use finishes...” (Mila i Canals et al. Mila i 
Canals 2007a). The modelling period is the time (years) over 
which the impacts caused by land use are integrated 
(i.e. difference between current impacts and a reference 
situation). 
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Considering an arbitrary and finite modelling period is a 
practical solution, which deviates from the recommendations 
in Mila i Canals et al. (2007a), who suggest that “the impacts 
on ecosystem quality should be assessed at least until a new 
steady state in ecosystem quality is reached by natural or 
human-induced relaxation” (regeneration); such modelling 
would introduce inconsistencies with different land transfor¬ 
mations, and would be hardly feasible in background systems 
(e.g. for LCI databases). The modelling period only has an 
effect on results when full recovery to reference situation is 
not achieved within this period (see Fig. 2). 

Taking a modelling period out of the three (20, 100 and 
500 years) applied by IPCC for C0 2 equivalents is a rea¬ 
sonable starting point. The modelling period should be 
chosen long enough to fully include the main transformation 
impacts measured by regeneration time. They can be esti¬ 
mated on the basis of observed re-vegetation on old 
quarries, past agricultural sites and closed railway tracks 
with known date of abandonment. As shown above, natural 
regeneration in temperate or boreal forest biomes takes 
clearly more time than 100 years. That is why a modelling 
period of 500 years was used by Muller-Wenk and Brandao 
(2010). In this guideline, a 500 years modelling period in 
land use LCIA is proposed as it is adequate in view of the 
long-term natural processes involved. 

3.3.2 Uncertainty 

Uncertainty is inherent in all steps of the calculation of global 
ecological impacts of land use. Firstly, using statistical mea¬ 
sures quantifying the statistical uncertainty of CFs is highly 
recommended. The first and best approach is to compare 
differences in the means of the CFs with statistical methods 
like ANOVA. The standard error (calculated as standard devi¬ 
ation divided by the root of sample size n) allows users to 
identify CFs with a high standard deviation in relation to the 
sample size, which are not appropriate for the use in LCA 
studies (Koellner and Scholz 2008). In cases where assump¬ 
tions of normal distributions cannot be met and parametric 
statistics can therefore not be applied, nonparametric statistics 
based on differences in median and percentiles should be used. 
Such information can be used in Monte Carlo simulations or 
sensitivity analysis to assess the overall uncertainty of the 
impact assessment. This helps users to interpret the current 
state of knowledge on land use impacts assessment (e.g. which 
impact pathways are most uncertain and which parameters are 
uncertain and need to be better investigated) and should also 
allow identifying research priorities. Second, in order to assess 
the data quality, we recommend to develop a pedigree matrix 
similar to what was presented in Ecoinvent 2.0 (Frischknecht 
and Jungbluth 2007, p. 45) and Ecoinvent 3.0 (Weidema et al. 
2011, p. 83) to assess the reliability of the data and the models 
behind specific characterisation factors, the completeness of 


the impact pathway as well as temporal/geographical 
representativeness. 

4 Discussion: problems of global impact assessment 

4.1 Simplifications and assumptions 

Assessing the ecological impacts of land use on a global 
scale is a complex matter. Several simplifications, assump¬ 
tions and decision choices are linked with this approach of 
the calculation of generic CFs for land use impacts. The 
main ones related to land transformation and occupation are 
the following: 

• Discrete land cover types based on environmental fac¬ 
tors like vegetation cover are the basis of assessing 
ecological impacts; however, estimates of ecosystem 
service changes may be inaccurate because species dis¬ 
tributions and their functional properties can vary con¬ 
siderably within one land cover type (Eigenbrod et al. 
2010). Replacing land cover maps with continuous en¬ 
vironmental information including species functional 
traits as input into ES models is likely to increase the 
model validity (see Lavorel et al. 2010). 

• Ecosystem quality Q remains constant over occupation 
time. This is not necessarily the case as during occupa¬ 
tion, Q can change (Mila i Canals 2007a). 

• Time and area of occupation can be fully substituted to 
achieve a constant output (e.g. one ton of crop). This 
might be true for the inventory, but in terms of ecolog¬ 
ical impact, it can matter either to use a small area for a 
long time or a large area for a short time (Koellner and 
Scholz 2007), or one single large plot versus many 
smaller ones. 

• Transformation time is zero. This assumption is made, 
because in most LCA applications, transformations im¬ 
ply worsening ecosystem quality, which in many cases 
requires not much time (e.g. deforestation and building 
of infrastructure). However, in cases of transformations 
improving the ecosystem quality (e.g. from agriculture 
to forest), this might not be appropriate. 

• Ecological impact is linearly increasing with the interven¬ 
tion of occupation and transformation. However, biodi¬ 
versity and ecosystem services might respond non- 
linearly to land use pressure (Carpenter et al. 2009). In 
LCA, a marginal approach is taken to account for this 
with the assumption that the “concentration” of land use 
is not considerably changed in the background system 
given the intervention at stake. However, the question is 
at which threshold does an intervention influence the 
start of background system and non-linear system 
response. 


<£) Springer 



Int J Life Cycle Assess (2013) 18:1188-1202 


1199 


• Biodiversity and the multiple ecosystem services are 
independent. However, research shows clearly there is 
an interaction between them (Balvanera et al. 2006; 
Bennett et al. 2009; Raudsepp-Hearne et al. 2010). A 
first step towards addressing impacts of biodiversity 
changes on ecosystem services is the assessment of 
functional diversity in LCA. It is also important to assess 
the relationship between impact results of different eco¬ 
system services (see e.g. Mila i Canals L et al. 2013). 

• Multiple drivers of ecosystem services and biodiversity 
loss do not interact. In most product systems, multiple 
environmental impacts occur, such as climate change, 
pollution, land use and the combined effects can lead to 
non-additive reactions of ecosystem services and 
biodiversity. 

The following assumptions are made with respect to 
ecosystem regeneration, which is a major factor to assess 
transformation impacts: 

• Regeneration is linear and independent from the 
land use history , i.e. only the last land use before 
abandonment is important and time of occupation is 
not relevant. This is certainly not appropriate for many 
cases as land use history determines the potential of 
ecosystems to recover. Ecosystems might also react 
non-linearly in their regeneration depending on their 
resilience (Folke et al. 2004). 

• Configuration of a landscape is not influencing regen¬ 
eration time and the degree of regeneration (e.g. perma¬ 
nent impacts). Factors like the proximity to primary 
habitat, connectivity of landscape and availability of 
remaining primary habitat are not considered for biodi¬ 
versity impacts. This could be addressed by a statistical 
approach as shown for Switzerland (Koellner 2003, p. 
97-125, Koellner et al. 2004). 

• No active restoration was assumed. Until now, we as¬ 
sumed that ecosystems regenerate naturally without hu¬ 
man influence. However, in some cases, technical means 
might contribute to ecosystem restoration (as for exam¬ 
ple applied in the restoration of mining sites) and can 
shorten the regeneration time. It is uncertain to which 
extent this is possible for ecosystem services and for 
biodiversity on larger scales, and particularly for land 
use in the background system. 

4.2 Separation of land use and water use 

The impact category land use is often closely linked to the 
water use category. This is typically the case for irrigated 
agriculture, where the irrigation is part of the land use 
practice, but at the same time linked to the water used. Land 
use changes can also induce indirect changes in the water 


cycle. For instance, deforestation or surface sealing can 
change water infiltration, evapotranspiration and runoff, 
leading to changes in groundwater aquifer recharge and 
thereby influence freshwater availability in a specific region. 
This eventually may cause adverse impacts on downstream 
ecosystems and their biodiversity (Pfister et al. 2009), which 
may have significant contributions to the overall system’s 
water-related impacts (Mila i Canals et al. 2010). More 
discussion on overlap and complementarities between the 
two impact categories is given in Koellner et al. (2011a). 

4.3 Support to decision-making 

It is acknowledged that the value of the LCA results for 
decision-making depends heavily on the appropriateness of 
the LCI and LCIA approaches taken. In this sense, it is 
important to highlight that the choice of reference situation 
is a key element for the decision maker and should be fully 
appreciated for a proper interpretation of the results. The 
approach suggested in this guidance is to consider a biome- 
dependent (quasi-)natural land cover as a reference. However, 
policy makers wanting to protect the current environmental 
quality may not be interested in knowing how far we are from 
an idealistic, quasi-natural situation. For this purpose, ap¬ 
proaches relying on realistic scenarios of ecosystem develop¬ 
ment on a landscape scale should be developed. In addition, 
the spatial disaggregation to biome level is probably relevant 
as a first step of spatial differentiation for global life cycles as 
shown in some case studies (see Mila i Canals L et al. 2013); 
however, if the decision maker wanted to support e.g. land use 
policy in a sub-biome level (e.g. in a country), then the LCI 
information (and LCIA CF) provided would need to be at a 
further level of detail. This also requires to clearly identify 
complementary methods at the landscape scale such as envi¬ 
ronmental impact assessment and regional ecosystem service 
models (e.g. Nelson et al. 2009). 

5 Conclusion and recommendations 

The principles and recommendations presented here are 
applied in the other papers of this special issue for the 
development of CFs for land use impacts on biodiversity 
and ecosystem services on a global scale. Such recommen¬ 
dations go beyond the general framework presented by Mila 
i Canals et al. (2007a) and explore further specific aspects 
such as allocation of transformation impacts over 20 years; 
pragmatic consideration of irreversible impacts with finite 
modelling period of 500 years; classification of impact 
pathways aligned with the Millennium Ecosystem Assess¬ 
ment (2005); discussion on the validity of inherent assump¬ 
tions in land use impact modelling (Section 4), consistent 
classification of different land use typologies and levels of 
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biogeographical differentiation (Koellner et al. 2013) that 
goes beyond the “archetypical land uses” suggested in Mila 
i Canals et al. (2007a). Table SI (Electronic Supplementary 
Material) summarises the impact pathways covered in the 
different papers and how the guideline presented in this 
paper has been followed in each one of them. 

The presented CFs provided by the other papers of this 
special issue are developed for the global assessment of direct 
land use impacts in the background system. For a more elab¬ 
orated assessment of land use impacts in the foreground sys¬ 
tem, regional models should be further developed and applied 
for biodiversity (Geyer et al. 2010a; Geyer et al. 2010b) and 
ecosystem services (Beck et al. 2010; Bare 2011). Indirect land 
use in a consequential approach is also not yet sufficiently 
considered in LCA, and even though it is part of the inventory, 
it may have important consequences in the impact assessment 
and needs further research. Brandao (2011) offers some key 
insights on the importance of indirect land use on the impact 
assessment results. 

Weighting of ecosystem services midpoints and their 
independence are not yet addressed. As they link to the 
human well-being according to the Millennium Ecosystem 
Assessment, such weighting can be based on the societal 
demand. For example, the relative importance can be 
expressed in their economic value (Costanza et al. 1997). 
Before such weighting is proposed, it would be worth to 
thoroughly explore whether the plethora of ecosystem ser¬ 
vices quantified in this special issue provide differentiated 
results, or whether they are all so inter-linked that assessing 
one of them would be enough. However, weightings of 
ecosystem services midpoints as well as many decisions 
done along the development of LCIA are clearly a matter 
of value choices. Further developments should make them 
consistently transparent for archetypical cultural perspec¬ 
tives (Hofstetter 1998; De Schryver et al. 2010). 
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